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ABSTRACT 

One of the most challenging issues for future wireless 
communication systems is the provision of quality-of- 
service (QoS) guarantees to users over the harsh 
wireless channels given the limited resource 
availability. The aim of this paper is to probe the 
problem of dynamic subcarrier and bit allocation in 
downlink of MIMO-OFDMA systems. Using Singular 
Value Decomposition, MIMO fading channel of each 
subcarrier is transformed into an equivalent bank of 
parallel SISO sub-channels. NSGA-II is handled 
which is a multi-objective Genetic Algorithm, for 
joint allocation of bits and subcarriers, in the 
downlink of MIMO-OFDMA system. NSGA-II is 
intended for optimization problems involving multiple 
conflicting objectives. Here the two conflicting 
objectives are Rate Maximization and Transmit Power 
Minimization. 

Keywords: MIMO-OFDM, Non-dominated Sorting 
genetic algorithm, Multi-objective optimization, 
Elitism, Singular value decomposition 

1. INTRODUCTION 

Wireless communications has emerged as one of the 
largest sectors of the telecommunications, industry, 
evolving from a niche business in the last decade to 
one of the most promising areas for growth in the 
recent era. 

The increasing demand for wireless multimedia 
services are in need of reliable and high-rate data 
communications over a wireless channel. High 
channel capacities or transmission rate for the same 
bandwidth with no additional power requirements by 
employing multiple antennas (MIMO) at the 
transmitter and receiver is contributed by Spatial 
multiplexing.(OFDM) efficiently uses the spectrum 


by spacing the channels closer, and it also has the 
ability to reduce 1ST The blend of these two 
technologies has been researched for the most 
hopeful applicant technique for the next generation 
wireless systems. Due to the frequency-selective 
fading, a major portion of the subcarriers assigned to 
a specific user would undergo severe fading and 
power will not be efficient to carry any data. As the 
number of user increases the Multiuser system 
creates channel. Adaptively assigning resources to 
each user provide a better performance compared to 
fixed schemes and it is based on their channel 
condition which can, abuse multiuser diversity [2], 

The MIMO/OFDM system enjoys more degrees of 
freedom for efficient resource allocation and it is due 
to the spatial parallelism and frequency selectivity of 
the channel. 

2. MIMO-OFDMA 

Multi-user MIMO-OFDMA is regarded as an 
important technology for increasing the flexibility and 
efficiency of wireless communication systems [2] .In 
the MIMO-OFDMA downlink system model from the 
base station (BS), all channel state information of 
each transmit and receive antennas is sent to the 
subcarrier and power allocation block through the 
feedback channels from all mobile users. The resource 
allocation scheme is forwarded to the MIMO-OFDM 
transmitter. The transmitter then loads each user’s 
data onto its allocated subcarriers. 

The resource allocation scheme is updated as soon as 
the channel information is composed and also the 
subcarrier and bit allocation information is guided to 
each user for the purpose of detection, through a 
separate channel [6], 
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MIMO exploits the space dimension to progress 
wireless systems capacity, range and consistency. It 
offers significant escalations in data throughput and 
link range without supplementary bandwidth or 
increased transmit power. MIMO achieves this goal 
by spreading the same total transmit power over the 
antennas to achieve an array gain that progresses the 
spectral efficiency (more bits per second per hertz of 
bandwidth). 

3. SYSTEM MODEL 

The system under deliberation is downlink of a 
MIMO-OFDMA system, with K users and N 
subcarriers. At the transceiver the perfect channel 
state information is assumed. Then each subcarrier 
composed of a narrowband channel with M antennas 
at both the receiver as well as the transmitter, which 
can be modeled by Mx M channel matrix and it is of 
H = [hi,j], where hi,j is the gain of the channel at the 
receive antenna i from transmit antenna j.The SVD of 
H is then written as: 

H=U A 1/2 V t (1) 

Where U and V are the unitary matrices, V’ 1 " denotes 
the transpose conjugate of V, and A is a diagonal 
matrix. The elements of A=diag([X L ....7 M ]) are real 
and they are arranged as ^^....^M^O.The SVD 
analysis is shadowed by the breakdown of channel 
matrix H into a number of self-governing orthogonal 
excitation modes, and it is known as eigen-modes of 
the channel. 


4. OPTIMIZATION PROBLEM 

The transmit power for the k t h user and the m t h Eigen¬ 
mode of the n* subcarrier must be equal to: 


Pk,n,m fk( c k,n,m)/ ^k ,n,r 


( 3 ) 


Where 


Ck,n,m 


is the bits allocated for unit channel 


gain. The overall transmit power is given by 
k n m f 

P n k ( C k,n,m) 

p T=2 J 2A»Z“ 

k=ln=l m=l 


Hc,n,m 


Where the subcarrier allocation indicator p kn is 
defined as 

Pk,n 

_ ( 1, if the nth subcarrier is allocated to kth user 
[0, otherwise 


Problem is framed to maximize the minimum data 
rate among entire users subject to the constraint that 
the total power cannot exceed a given value. The 
power objective is modified little and assumes that 
the total available transmission power is limited to a 
certain range with a typical value P|.Second 
objective to bring the total power as close as possible. 
The multi-objective optimization problem is: 

N M 


max min R k = max min 

c k,n,m 'Pk,n k c k,n,m,p k n k 


min 

Ck.n.m ,p kn 


II c k,n,m Pk,n 

n- 
M 


71 = 1 772=1 

* * M f k (<*,„.„) 


/C — 1 72=1 


772 = 1 


l k,n,m 


The transmit preceding and receiver shaping convert 
the MIMO channel into R H (R H is the rank of the 
matrix H) parallel SISO channels with input x and 
output y 

Since from the SVD,it has, 
ym V^ m n m> form—1, M (2) 


Subject to constraints: 

K 


NMD r 


K =1 
K 

C 2:^P k ,n = 1 ’ andp kn = {0,l} 


K=1 


Here n = U^n ,and n denotes to Additive White 
Gaussian Noise(AWGN). 

In the frequency selective fading different subcarrier 
constitute of different MIMO channels. The MIMO 
channel matrix H kn for the /causer and the 
n th subcarrier is disintegrated into M parallel SISO 
channels, where channel gains are given by 

V^k,n,m(^k,n,l-- • • --^k,n,M-0)- 


Rk denotes the transmit bits per OFDM symbol for 
the k t h user, D is a set of available bits in the adaptive 
modulation, D max denotes the maximum value in the 
set D. Constraint C 2 takes care of the fact that one 
user is assigned with one subcarrier only. 

5. NSGA-II 

Non-dominated Sorting Genetic Algorithm is to use to 
progress the adaptive fit of a population of candidate 
solutions to a Pareto front which is inhibited by a set 
of objective functions [9], The algorithm uses an 
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evolutionary process with for evolutionary operators 
that include selection, genetic crossover, and genetic 
mutation. The population is sorted into a hierarchy of 
sub-populations based on the ordering of Pareto 
dominance. Similarity between members of each sub¬ 
group is evaluated on the Pareto front, and the 
resulting groups and similarity measures are used to 
endorse a various front of non-dominated solutions. 
NSGA-II algorithm includes the following steps to be 
employed: 

5.1 Population initialization 

The number of individuals in the population P and the 
number of generations G are initialized and can be 
varied for different runs of the algorithm. Each 
individual is created by generating 0’s and l’s 
randomly. 

5.2 Evaluate the objective functions 

For each individual the fitness values of the objective 
function are calculated. 

5.3 Non-dominated sorting 

Based on non-domination the population is sorted into 
fronts and each individual is assigned a rank. After the 
completion of non-dominated sorting the crowding 
distance is assigned.The non-dominated sorting 
algorithm is done by following steps: 

Step 1 : From the main population P each individual p 
is made to do the following: 

• Initialise Sp= <j> 

{set of individual being dominated by p} 

• Initialise Np=0{ no of individual that dominates 
Pi- 

• Consider individual p in population P 

> If p dominates q then S p is 
S p ={qUSp} 

> Else Np=Np+l 

Step 2: IfNp=0 (no individual dominates p and it 
belongs to first front).Set p ran k=l.First front update is 
made with addition of p to the set(Fi={ p U Fi}). 

Step 3: Intialize front counter i=l 

(Carry out for the entire individual in main population 

P) 

• For each individual q in Sp 

> Nq=Nq-l (Nq is domination count) 

> If Nq=0 then 

> q rank = i+l 

• update Q ={ Q U q } 


Step 4\ If Fi not equal to 0 then 

Q=0{set for storing individual in i+l front} 

For each individual p in Fi 

Step 5\ Increment the front counter Fi = Fi + 1 

Step 6:Set Q in the next front Fi=Q. 

5.4 Crowding distance calculation 

In addition to fitness value for each individual 
crowding distance is calculated. The crowding 
distance is calculated as: 

Step 1: For each front Fi, in a non-dominated set 1,1 is 
the number of solution 

Step 2: For all the individuals, initialize the distance to 
0 i.e. ,Fi(dj) = 0, j corresponds to j t h individual in 
Fi.And for each objective function m 

Step 3: Based on the objective m sort the individuals 
in front Fi (I=sort(Fi,m)). 

Step 4 : Boundary solution is assigned to infinite 
distance and for each individual in front Fi(I(dl) = oo 
and I(dl) = oo).By this for all other points the 
boundary points are always selected. For j = 2 to (1-1) 

Step 5: For j = 2 to (1-1) 

• I(dj) = I(dj) + (I(j+ l).m - I(j-l.m)/fmax(m) - 
fmin(m) 

• I(j)m is the value of the jth individual of the mth 
objective function in I. f max(m) and f min(m) are 
the maximum and minimum values of the mth 
objective function. 

Based on the rank and crowding distance parents are 
selected from the population by using binary 
tournament selection. 

5.5 Crossover and Mutation 

Crossover is the mechanism by which design 
characteristics between any paired individuals are 
exchanged to form two new (child) individuals and 
single point crossover the data beyond the point in 
either organism string is swapped between the two 
parent organisms. 

Mutation is a genetic operator used to maintain 
genetic diversity from one generation of a population 
of chromosomes to the next. The probability of 
mutation is given as P m and it is usually 100 times less 
than P c , where P c is the probability of crossover. 
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5.6 Generation of New Population 

The new population is filled up by taking the best 

individuals from the combined population based on 

rank and crowding distance. 

6. SIMULATION RESULTS 

Simulation results for permutations of condition in 

Table I. 

A. Plot of Power optimization function that minimize 
the difference between the typical Power available 
(PT) and the Power consumed (Ptotal) by the 
current allocation scheme for 2 users 

B. Plot of rate optimization function which is 
1/(Minimum rate + 1) for 2 users 

C. Plot of Power optimization function which is 
minimizing the difference between the typical 
Power available (PT) and the Power consumed 


(Ptotal) by the current allocation scheme for 4 
users 

D. Plot of rate optimization function which is 1/ 
(minimum rate + 1), for 4 users 

E. Plot of Power optimization function which is 
minimizing the difference between the typical 
Power available (PT) and the Power consumed 
(Ptotal) by the current allocation scheme for 8 
users 

F. Plot of rate optimization function which is 
1/(minimum rate +1), for 8 users 

G. Plot of Power optimization function which is 
minimizing the difference between the typical 
Power available (PT) and the Power consumed (P 
total) by the current allocation scheme for 16 users 

H. Plot of rate optimization function which is 
1 /(minimum rate + 1) for 16 users 

.El 


Maximum number of generations required for convergence for different set of experiments 


s . 

No 

Users 

(K) 

Subcarriers 

(N) 

Population 

Maximum 

generations 

Pc 

Pm 

Dmax 

PT 

(dB) 

Maximum 

generation 

for 

convergence 

1 

2 

64 

50 

80 

0.8 

0.02 

8 

13 

35 

2 

4 

64 

30 

100 

0.8 

0.02 

16 

24 

28 

3 

8 

64 

50 

50 

0.9 

0.03 

8 

35 

26 

4 

16 

64 

30 

50 

0.9 

0.03 

8 

56 

25 



(a)Number of Generations vs abs(PT-Ptotal) 



(b) Number of Generations vs l/(min user 
rate+l)(bps) 



(c) Number of Generations vs abs(PT-Ptotal) 



rate+l)(bps) 
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rate+l)(bps) 




N=64,K=16 

Pc=0.9,Pm=0.03 

PT=56dB,Dmax=8 



3 



Number of Generations 


(g) Number of Generations vs abs(PT-Ptotal) 



(h) Number of Generations vs l/(min user 
rate+l)(bps) 

Figl. Simulation results for permutations of 
conditions given in Table I 



Fig 2 sum rate (b/s hz) vs number of users for 
snr=24.6db 



Fig3. Sum rate (b/s hz) vs number of users for 
snr=4.6db 


The Simulation results from FIG 2 and FIG 3 indicate 
the optimized Sum rates (b/s/Hz)for fixed number of 
sub-carriers(N=64)versus number of users(K=4,8,16) 
for SNR=24.6Db (Dmax=16 bits) (Fig 2) and for 
SNR=4.6db (D max =8bits) (Fig 3). 

7. CONCLUSION 

NSGA-II has been proposed for subcarrier allocation 
in MIMO-OFDM system for the multiuser and 
simultaneous bit-loading .The solution generated by 
the algorithm is found to be suitable for different sets 
of users and subcarriers taking multiple conflicting 
objectives into the algorithm and it guarantees 
fairness among the users by means of optimizing the 
rate. 
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